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Abstract 

Super Cooled Large Droplet (SLD) collection efficiency calculations were performed for several 
geometries using the LEW1CE3D Version 3 software. The computations were performed using the NASA 
Glenn SLD splashing model which has been incorporated into the LEW1CE3D Version 3 software. 
Comparisons to experiment were made where available. The geometries included two straight wings, a 
swept 64A008 wing tip, two high lift geometries, and the generic commercial transport DLR-F4 wing 
body configuration. In general the LEW1CE3D Version 3 computations compared well with the 2D 
LEW1CE 3.2.2 results and with experimental data where available. 

Nomenclature 


AOA angle-of-attack 

AAOA aircraft angle-of-attack 

C p pressure coefficient 

c chord 

d droplet diameter, pm 

K l LEW1CE splashing parameter 

LWC Liquid Water Content, g/m 3 

MVD Median Volume Diameter 


Introduction 

Over the last several years work has been ongoing to develop tools to analyze aircraft configurations 
subjected to Super-Cooled Large Droplet (SLD) conditions (Ref. 1). These conditions, which are outside 
the FAA Appendix C (Ref. 2) certification envelope have been responsible for several aviation accidents. 
New certification rules are pending which will require aircraft to fly safely through these conditions. 
These new rules will generate new requirements for aircraft ice protection system design and certification 
and for the tools which aid in these processes. 

One of the major challenges in developing tools for SLD condition has been the quantification of the 
effect of splashing of large drops. The larger drops contained in the SLD environment (> 50 pm) carry 
more momentum and arrive at higher speeds at the aircraft surface then the smaller drops defined within 
the Appendix C envelope. Surface tension for the large SLD drops is not sufficient to hold the drops 
together under the high impact forces as it is for the smaller drops (« 20 pm), which results in droplet 
breakup and splash. Droplet breakup up and splash can result in a substantial amount of mass being 
ejected from the surface. 

Ignoring the effects of splashing for SLD conditions in a computational simulation can result in an 
over-prediction of collection efficiency. This over-prediction, although conservative, can in many cases 
yield intractable designs. For this reason, tools which accurately predict the effect of splashing are needed 
to aid in the design, development and certification of aircraft systems for SLD conditions. 
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NASA Glenn has incoiporated SLD splashing models into several of their icing tools to aid in the 
design, development, and certification of aircraft ice protection systems. These tools include the two- 
dimensional (2D) LEW1CE 3.2.2 software (Ref. 3), hereafter referred to as LEW1CE, and the three- 
dimensional (3D) LEW1CE3D V3 software (Ref. 4), hereafter referred to as LEW1CE3D. The LEW1CE 
software, which contains the original NASA SLD splashing model, has been previously validated against 
large droplet collection efficiency data gathered during NASA-Wichita State University (WSU) Icing 
Research Tunnel (1RT) droplet impingement tests (Ref. 5). The work presented in this paper represents 
the validation of LEWICE3D with the newly incoiporated NASA Glenn SLD model. 

Collection efficiency results are presented for two 2D wings, a swept 64A008 wing tip, two high lift 
models and a generic commercial transport. The wings include an MS-317 (Ref. 6) and an NLF(1)-0414 
(Ref. 7). The high lift models include a straight, three element high lift wing model developed by 
McDonnell Douglas Aerospace (Ref. 8), and the generic trapezoidal high-lift wing (trap wing) developed 
by NASA (Ref. 9). The DLR-F4 wing body model developed as a CFD test bed was used to represent a 
commercial transport configuration (Ref. 10). Experimental SLD data was available for the 2D wings and 
for the three element wing model. No SLD data exists for the trap wing and for the DLR-F4 
configuration. 


Experiment 

Experimental Apparatus 

The aerodynamic and collection efficiency tests were carried out in the NASA Glenn 1RT. The test 
equipment included the 1RT (Fig. 1), the ESCORT data system, a special spray system for the 
impingement tests (Fig. 2), a laser reflectometer for impingement efficiency data reduction (Fig. 3), 
blotter strips attached to the models (Fig. 4), and the four test models (Figs. 5(a) to 8(a)). 

The IRT facility can provide a range of airspeeds, angles-of-attack, temperatures, liquid water 
contents (LWC), and drop sizes (Ref. 1 1). The IRT has a 2.47- by 1.82-rn test section with a maximum 
airspeed of 134 m/s (empty tunnel). Angle-of-attack is controlled by a movable turntable to which the 
models are mounted. A refrigeration system allows year-round testing at temperatures from -29 to 10 °C. 
The spray system located upstream of the test section can provide a cloud with a range of LWC of 0.25 to 
3.0 gm/m 3 and a median volumetric droplet diameter (MVD) size range of 12 to 200 pm. 

The ESCORT system was developed at Glenn to aid in storage, processing, and analysis of large 
amounts of data (e.g. temperature, pressure) produced in various experiments at Glenn Research Center. 

In this test ESCORT was used to store tunnel total temperature, total pressure, free stream airspeed, 
surface pressure, produce real time calculations and display pertinent parameters. The storage sequence 
for each data point was initiated by the researcher in the control room. A separate program was used to do 
a more complete post run analysis. 

The spray requirements for the impingement tests precipitated the need for a different spray system 
(Fig. 2) than is available in the IRT (Ref. 5). The IRT spray system could not produce the short (2 to 
9 sec), stable sprays (i.e. constant LWC and drop size) required to prevent blotter strip saturation. There 
were also concerns that the dye would contaminate the IRT spray system. The new system consisted of 1 6 
nozzles and a supply tank located at the IRT spray bar station (Fig. 2). The system featured short supply 
lines which enabled short, stable sprays. 

One unique feature of the current technique is the laser reflectometer used to determine the local 
collection efficiency (Fig. 3). The device measures the local reflectance of the blotter strip and correlates 
this to the local collection efficiency (Ref. 5). The device saved considerable time in the data reduction of 
the blotter strips. 

The MS-3 17 airfoil was designed in the mid 1970s to bridge the gap between low speed and 
supercritical airfoils for application on general aviation aircraft. The airfoil was designed for a lift 
coefficient of 0.3, Reynolds number of 14 million and a Mach number of 0.68. The IRT MS-3 17 model 
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(Fig. 5(a)) was an unswept, untapered, full span model. The chord was 0.9144 m, and the span was 
1.8288 m. The model contained 49 pressure taps located along the wing centerline. 

The natural-laminar-flow NLF(1)-0414 wing model was an unswept, untapered, full span model built 
for the 1RT (Fig. 6(a)). The NLF(1)-0414 airfoil was designed in the early 1980s as a medium speed 
airfoil with low section drag and high maximum lift. The airfoil was designed for a Reynolds number of 
10 million, a Mach number of 0.4, and a lift coefficient of 0.4 with 0.7c laminar flow on both the upper 
and lower surfaces. The wing had a chord of 1.2192 m in. and a span of 1.8288 m. The pressure 
instrumentation included 124 taps located along the center span of the model. 

The swept, tapered NACA 64A008 wing tip represents a typical business jet horizontal tail 
(Fig. 7(a)). The model had a span of 1.2192 m, a leading edge sweep angle of 29.1°, a root chord of 
1.1621 m, a taper ratio of 0.62 and a cylindrical tip. The model had two rows of 31 pressure taps at the 
50 and 90 percent span stations. 

The multi-element 1RT model was full span (1.8288 m) and had a nested chord length of 0.9144 m 
(Fig. 8(a)). The aluminum model had three elements (slat, main element, and flap) and contained 128 
pressure taps. The taps were located along the model centerline and along two span-wise rows on the 
upper surface and on one span-wise row on the lower surface. 

Experimental Testing 

Two types of testing were done in the 1RT: aero-performance and impingement efficiency testing. 

The aero-performance testing involved the use of the NASA Glenn ESCORT system to take surface pres- 
sure measurements. The impingement efficiency testing involved the use of a dye tracer technique to 
measure the location and amount of water impinging on the model. 

Surface pressures were measured on the wing models using the ESCORT system. Pressure measure- 
ments were taken at an airspeed of 75 m/s, at angles-of-attack of 0° and 8° for the MS-3 17 model, 0° and 
8° for the NLF(1)-0414 model, 0° and 6° for the swept 64A008 wing model, and 0° for multi-element 
model. These pressure measurements were used to generate coefficient of pressure distributions for the 
comparisons presented in this paper. 

The experimental technique used in the current tests to determine the impingement characteristics of a 
body is one that was developed in the early 1950s with a few modifications (Ref. 5). The technique 
involved spraying a dye -water solution of a known concentration onto a model covered with blotter strips. 
Figure 4 shows a typical blotter installation for the swept 64A008 wing. The result was that the local 
impingement efficiency on the blotter strips was related to the variation in color intensity. That is, the 
areas of higher impingement rate are darker and those with lower impingement rate are lighter. 

Several steps were necessary to prepare the 1RT for impingement testing. The specially designed 
spray system had to be installed and adjusted to produce a uniform cloud. The local LWC had to be 
measured at each blotter strip location (with the tunnel empty) for every spray and tunnel condition to 
account for any cloud non uniformity that existed after the final spray adjustment. After these adjustments 
and measurements were made the model was inserted and tested. Each test point was repeated four times 
to obtain a measure of the repeatability (Ref. 8). 

A typical test point involved several steps. The model was cleaned and blotter strips were attached at 
points of interest. The spray was then made, the blotter strips were removed, and labeled, and the model 
was cleaned and made ready for the next condition. 

The models were tested for several MVD drop sizes ( 1 1 -, 2 1 -, and 92-pm) and angles-of-attack. The 
work in this paper will focus on the 92 pm MVD drop size case which is indicative of an SLD condition. 
Collection efficiency data was taken at an airspeed of 80 m/s, a static temperature of 277 K, and static 
pressure of 94,540 Pa. The angle-of attacks tested were 0° and 8° for the MS-3 17 model, 0° and 8° for the 
NLF(1)-0414 model, 0° and 6° for the swept 64A008 wing model, and 0° for the multi-element model. 
The drop size tests were repeated four times for each condition to generate a statistical average. 

Measurements of the spray cloud drop distributions were made using several drop size measuring 
probes. The NASA Glenn Forward Scattering Spectrometer Probe (FSSP) and the Optical Array Probe 
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(OAP) were used to measure the droplet concentrations and MVD’s for the various spray conditions. 
These measurements were used to generate the 10 bin and 27 bin distributions used for the analysis. 
These distributions are shown in Tables 1 and 2. 


Analytical Method 

A variety of 2D and 3D flow and icing codes were used for the analysis. Some of the flow solutions 
were generated at NASA Glenn while others were provided by industry. The LEW1CE and LEW1CE3D 
programs were used to generate the collection efficiency results. 

CFD++ 

The CFD++ flow solver was developed at Metacomp Technologies (Refs. 12 and 13). CFD++ is a 
finite volume, unsteady, Navier-Stokes flow solver which can accommodate unstructured, hybrid, and 
overset grids and can generate flow solutions for flow regimes covering incompressible to supersonic. 

The solver can handle a variety of flow problems such as high speed reacting/non reacting flows, multi 
phase flow, conjugate heat transfer and electrostatics. The software supports a number of popular 
turbulence models and can be run on a variety of parallel platforms including, Linux, Unix and Windows. 

PMARC 

PMARC is a first order 3D potential flow panel code (Ref. 14). Geometries are represented as 
quadrilaterals which have constant doublet and source distribution. The formulation used in PMARC 
results in a solution that is second order accurate allowing for accurate flow solutions with fewer panels 
and less CPU time than other first order methods. The disadvantage of this method is that because a 
numerical differentiation is used to generate the velocity distribution careful panel modeling is required to 
prevent numerical errors. The code can generate solutions for internal and external compressible flows 
and can handle very large problems (-10,000 panels). 

USM3D 

USM3D is a 3D compressible, unstructured, Navier-Stokes flow solver (Ref. 15). The solver 
developed by Frink is a tetrahedral, cell centered, finite volume based flow solver. The parallel solver 
which can solve steady or unsteady cases employs local time stepping and Residual smoothing to aid in 
convergence. 


NSU3D 

NUS3D is a 3D unstructured, compressible, Navier-Stokes flow solver (Ref. 16). The parallel solver 
developed by Mavriplis employs a Roe-Rieman solver and supports multi-grid and low Mach number 
pre-conditioning to aid convergence. 


INS2D 

1NS2D is a 2D, multi-block, structured Navier-Stokes flow solver which supports overset grids 
(Ref. 17). The code developed by Rogers can handle steady and unsteady problems. 

LEWICE 3.2.2 

LEW1CE 3.2.2 is a 2D, multi-time step ice accretion code (Ref. 3). The code uses a 2D panel code to 
generate a new flow-field at each time step. The latest version of the software enables the analysis of SLD 
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icing conditions and of hot air and electro-thermal anti-icing systems. LEW1CE performs its analysis for 
2D geometries in minutes on a desktop PC, allowing the user to run several parameter studies for design 
purposes. The ice shape predictions have been validated against a wide variety of experimental 
conditions. 


LEWICE3D Version 3 

The LEW1CE3D grid based code incorporates droplet trajectory, heat transfer and ice shape 
calculation into a single computer program (Ref. 4). The code can handle generic multi-block structured 
grid based flow solutions, unstructured grid based flow solutions, simple Cartesian grids with surface 
patches, and adaptive grids with surface patches. The latter two methods allow the use of generic panel 
code input which, when combined with LEW1CE3D, is a computationally efficient method for generating 
ice shapes. The code can handle overlapping and internal grids and can handle multiple planes of 
symmetry. Calculations of arbitrary streamlines and trajectories are possible. The code has the capability 
to calculate tangent trajectories and impingement efficiencies for single droplets or droplet distributions 
using a standard area-based collection efficiency method or a Monte Carlo based method. Ice accretions 
can be calculated at arbitrary regions of interest in either a surface normal or tangent droplet trajectory 
direction. The program can run on a variety of single processor and parallel computers, including Unix, 
Linux and Windows. 

Version 3 of the LEW1CE3D software, which incorporates several new features, was used for the 
analysis. These features include a new particle splash and bounce algorithm, a new geometry handling 
scheme which allows complex mirroring, transformation and relative motion of input grid blocks and a 
new algorithm which calculates block to block collection efficiencies. These new additions will enable 
users to analyze SLD conditions and to calculate collection efficiency with particle splash and bounce 
effects through turbo-machinery. 


Configuration 

The six models chosen for the analysis encompassed a variety of models and flow solution strategies. 
The LEW1CE code, which has its own internal flow solver, uses a 2D Hess-Smith based panel code. The 
LEW1CE3D code, which requires a grid based flow solution, employed solutions from the USM3D, 
1NS2D, CFD++ and PMARC. The 3D solutions involved various approximations and boundary 
conditions. 

The panel model for the LEW1CE MS-3 17 model contained 108 elements while that for the LEW1CE 
NLF(1)-014 model contained 105 elements. 

The USM3D flow solver was used to generate flow solutions for the MS-317 and NLF(1)-0414 
LEW1CE3D analysis. The flow solutions were inviscid and mirror boundary conditions were used at the 
wing tips. The mirrored wing tip boundary condition was used to generate an extruded, 2D flow-field 
approximation. The grids for the MS-3 17 and NLF(1)-0414 models were generated using the GridTool 
and VGR1D programs (Refs. 18 and 19) distributed in the Tetruss Package. The MS-317 grid model 
contained 5,153,753 tetrahedrons and 909,708 grid points. The NLF(l) grid contained 7,264,130 
tetrahedron and 1,273,226 grid points. The surfaces for the MS-317 and NLF(1)-0414 models are shown 
in Figures 5(b) and 6(b) respectively. 

The swept 64A008 wing flow model was generated using the PMARC/1CEGR1D panel code interface 
for LEW1CE3D. The 1CEGR1D program (Ref. 20) was used to generate a 3D “oct-tree” type grid about 
the wing and the PMARC code was used to calculate velocities at the grid nodes. PMARC was executed 
using a steady, isolated flow with a y = 0 plane of symmetry (Fig. 7(b)). The swept 64A008 panel model 
contained 2,303 panels and the grid contained 216,708 nodes. 

Two flow solver methodologies were employed for the multi-element wing. The 1NS2D program was 
used to generate an isolated 2D incompressible, viscous flow solution for the wing. The 2D, structured, 
multi-block, chimera based solution was extruded to generate a 3D solution. The 1NS2D extruded, 3D 
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grid contained four blocks and 725,706 grid points. The USM3D flow solver was also used to generate 
3D flow solutions for the analysis. The USM3D flow solutions were viscous, compressible and the flow 
model included the 1RT test section walls. Separate grid models were developed for the 0° and -1.5° 
angle-of-attack solutions. The 0° angle-of-attack grid model contained 53,268,615 tetrahedrons and 
9,164,492 grid points. The -1.5° angle-of-attack grid model contained 52,814,552 tetrahedrons and 
9,084,829 grid points. The surface model for the 1NS2D model is shown in Figure 8(b) while that for the 
USM3D model is shown in Figure 8(c). 

The CFD++ unstructured flow solver was used to generate the flow solution for the trap wing. The 
CFD++ solution was generated at a Reynolds Number of 8.8 million, a Mach number of 0.23, and an 
aircraft-angle-of-attack of 8° on a model which had a span of 6.27 m. The airspeed for this case was 
70 m/sec. The half-plane grid contained 3,804,849 nodes, 9,71 1,808 volume elements and 173,287 
surface elements. The volume grid was developed using the AFLR3 grid generation program (Ref. 21). 
The surface model for the trap wing is shown in Figure 9. 

The DLR-F4 flow solution was generated using the NSU3D flow solver. The half-plane flow solution 
was inviscid, and incompressible and was generated for an aircraft-angle-of-attack of 0°. The volume grid 
contained 1,305,296 volume elements and 257,240 nodes. The surface model for the trap wing is shown 
in Figure 10. 


Analysis 

Results from the comparison of measured and predicted pressure distribution and collection 
efficiency for the models tested in the 1RT will be discussed. Results for the trap wing and DLR-F4 
Airplane analysis, for which no experimental collection efficiency data exists, will also be discussed. 

The coefficient of pressure distribution comparisons for the MS-3 1 7 model are shown if Figure 1 1 . 
The measured values are depicted with the analytical results from the LEW1CE and USM3D flow 
solutions. The angle-of-attack used for both the LEW1CE and USM3D flow solutions was that which 
gave the best comparison to experiment. This was done to reduce the influence of flow-field deficiencies 
on the calculated collection efficiency results. Deficiencies in the inviscid, 2D nature of the LEW1CE and 
USM3D flow-fields and in flow angularity in the 1RT can result in as much as a 2° difference between the 
measured turn table angle and the locally predicted angle-of-attack. A 2° difference in angle-of-attack can 
cause a significant change in collection efficiency. Because the main intent of this work was to validate 
the calculation of collection efficiency and not the calculation of flow the adjustment of AOA to generate 
the best pressure distribution comparisons was considered a reasonable approach. For the 0° experimental 
condition the best comparison was at angle-of-attack of -2.0° and -1.85° respectively for LEW1CE and 
for USM3D. The angle-of-attack for the best comparison for the 8° experimental data was 6. 15° for both 
LEW1CE and USM3D. 

The calculated pressure distributions compare well to each other and to the experimental results for 
both the 0° and 8° angles-of-attack for the MS-3 17 model. The maximum coefficient of pressure is 
noticeably higher for the USM3D prediction than for the experimental and LEW1CE results for the 8° 
case. This is probably because the USM3D flow model had higher resolution near the pressure peak and 
was therefore able to pick up the pressure peak more accurately. The trailing edge coefficient of pressures 
are under-predicted for both the 0° and 8° cases. This is probably due to viscous effects such as a 
“thickening” of the trailing edge due to boundary layer growth or trailing edge separation which the 
inviscid flow solutions do not model. 

Collection efficiency results for the MS-317 model are shown in Figures 12 and 13. Results are 
shown with and without splashing and for both LEW1CE and LEW1CE3D. From the figures one can see 
that the collection models for the 3D and 2D codes agree well. The results for the 27 bin distribution and 
the 10 bin distribution also compare well with the 10 bin distribution for both the 0° and 8° angle-of- 
attack cases verifying that the 10 bin approximation for the SLD distribution is reasonable. Both 
programs over-predict the peak collection efficiency by about 18 percent for the 0° and 8° angle-of-attack 
cases and over-predict the point at which the collection efficiency goes to zero on both the upper and 
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lower surfaces (i.e. impingement limits). One can also see that SLD calculations agree well with the 
experimental data and that neglecting splashing results in a large over-prediction of the collection 
efficiency near the impingement limits. The splashing model results in an over-prediction of the 
impingement limits and an over-prediction in the collection efficiency near the mid span of both the upper 
and lower impingement regions. The over-prediction is more severe for the higher AOA case and 
approaches about 30 percent for the lower surface of the model. This over-prediction though only results 
in about a 5 to 10 percent over-prediction in total collection efficiency (area under the collection 
efficiency curve) which is much better than the non-splashing result for which the over-prediction is 1 5 to 
20 percent. The splashing result is conservative in that it is an over-prediction but is not overly 
conservative as is the non-splashing results. One can also see that the LEW1CE3D predicted results are 
slightly larger than the LEW1CE results near the impingement limits. This is because a splashing limiter 
was required by LEW1CE3D to keep trajectory integration times reasonable for some particles impacting 
at shallow angles. This limiter resulted in larger collection efficiencies for LEW1CE3D than for LEW1CE, 
which did not require a limiter, near the impingement limits where the impact angles are typically small. 

The LEW1CE3D splash limiter employed sets the minimum impact angle for a given LEW1CE 
splashing parameter, K L , (Ref. 22) for which a particle splash can occur. For a small impact angle, which 
is the angle between the surface and the impacting droplet, the NASA Glenn SLD model yields a small 
splash angle, which is the angle at which the splashed droplet rebounds from the surface. Small splash 
angles yield particle positions very close to the body and generate splash velocities nearly parallel to the 
surface. For large or waiped surface cells numerical issues in defining the exact location of the 3D surface 
can yield erroneous impacts and problematic velocity searches for particles which are very close to the 
surface. The result being that that the rebounding particles impact within a few integration steps of 
rebound. This process repeats itself with the particle skipping many times until it has exited the 
computational domain. Particle “skipping” dramatically increases trajectory integration time and, in 
general, results in no additional impingement. By employing a splashing limiter the droplet trajectories 
are terminated at the impact point for shallow impact angles which avoids particle “skipping” thus 
decreasing particle integration time. This limiter results in more collection at the impact location then if 
the particle were allowed to splash. 

Results for the NLF(1)-0414 model are shown in Figures 14 to 16. The coefficient of pressure 
comparisons for the 0° and 8° angles-of-attack are shown in Figure 14. As for the MS-317 case the 
analytical angles-of-attack used were those that gave the best comparison to the experiment pressure 
distributions. For the 0° experimental condition the best comparison was at angle-of-attack of -0.75° for 
both LEW1CE and for USM3D. The angle-of-attack for the best comparison for the 8° experimental data 
was 6.3° for both LEW1CE and USM3D. The analytical pressure coefficients compare well with 
experiment for the both the 0° and 8° cases. The pressure peak predicted is slightly higher than that 
measured and predicted by LEW1CE. This is attributed to the more highly refined surface and to the 
inviscid nature of the USM3D solution. As for the MS-3 17 model the trailing edge coefficient of 
pressures are under-predicted for both the 0° and 8° cases. This again is probably due to viscous effects 
which the inviscid flow solutions do not model. 

Collection efficiency results are shown in Figures 15 and 16 with and without splashing and for the 
both LEW1CE and LEW1CE3D. From the figures one can see that the collection efficiency models for the 
3D and 2D codes agree well. The LEW1CE3D non-splashing collection efficiency result shows a region 
of increased collection near the trailing edge due to the cusped nature of the trailing edge for the 8° angle- 
of-attack case. The LEW1CE and experimental results do not show this region of increased beta probably 
due to a small difference in the local angle-of-attack between the wind tunnel model and LEW1CE 
models. The results for the LEW1CE 27 bin and 10 bin calculations agree for both angles-of-attack. One 
can also see, as was the case for the MS-3 1 7 model, that the SLD calculations agree much better with the 
experimental data and that neglecting splashing results in a large over-prediction of the collection 
efficiency near the impingement limits. The collection efficiencies are slightly over-predicted from the 
mid-span to the trailing edge for both the 0° and 8° cases by LEW1CE3D due to the impact angle limiter 
mentioned above. As in the MS-3 17 case, the splashing model over-predicts the total collection efficiency 
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(< 1 0 percent) but the over-prediction is much less than that of the standard no-splash model 
(>20 percent). The splashing model results are conservative but not overly so as was the case for 
the non-splashing result. 

The results for the NACA 64-A008 wing are shown in Figures 17 to 19. Because the model was a 3D 
swept, tapered wing results are only presented for USM3D and LEW1CE3D. Angles-of-attack of 0° and 
6.3° were used to best simulate the flow for the 0° and 6° experimental cases respectively. 

The pressure distribution comparisons for the NACA 64A008 wing are shown in Figure 17. The 
USM3D results compare well with the experimental data except near the pressure peak and at the trailing 
edge. As was the case for the previous models, the higher resolution of the surface model results in ability 
to better resolve the pressure peak although this peak is probably over-predicted due to the inviscid nature 
of the flow solution. The trailing edge coefficient of pressures were under-predicted due to viscous effects 
not being captured by the inviscid flow model. 

The collection efficiency results for the NACA 64A008 wing are shown if Figures 1 8 and 1 9. The 
results for the standard model show a large over-prediction in local and total collection efficiency 
(> 20 percent). The collection efficiency results for the 10 bin calculation with splashing agrees well with 
the experimental data (< 5 percent). The splashing model over-predicts the peak collection efficiency by 
about 13 percent for the 0° and 6° angle-of-attack cases. For the 6° case there is an under-prediction of the 
collection efficiency near the mid-span on the pressure side of the wing and an over-prediction of the 
collection efficiency near the impingement limits. The impingement limits are also over-predicted for the 
6° angle-of-attack case. The splashing again produces a reasonably conservative result for the overall 
collected mass. 

Pressure and collection efficiency distributions are shown for the multi-element model in Figures 2 1 
to 26. Experimental pressure distributions are shown along with analytical pressure distributions for the 
1NS2D and USM3D based flow solution for several angles-of-attack. The 1NS2D solution was generated 
for an angle-of-attack of 0°. The USM3D solutions were generated for 0° and -1.5°. The -1.5° case was 
generated to better match the pressure distribution on the slat. 

Figures 21 to 23 display the pressure distribution comparisons for the multi-element model. The 
1NS2D results and USM3D computations for 0° angle-of-attack agree well except near the pressure peak 
for the main element and flap. The agreement of the USM3D and 1NS2D solutions at an angle-of-attack 
of 0° was considered quite remarkable considering that the 1NS2D solution was incompressible and that it 
did not include wall effects. The larger peak value for the main element and flap for the 1NS2D solution is 
partly the result of the incompressible assumption for the flow and due to the higher resolution of the 
1NS2D surface grid. The 1NS2D and USM3D pressure distribution for 0° angle-of-attack agree well with 
the experiment for the main element and flap but the comparison for the slat was considered poor. The 
comparisons imply an over-prediction in local angle-of-attack for the slat and an under-prediction in 
angle-of-attack for the flap. An attempt was made to improve the agreement by running a USM3D flow 
solution at an angle-of-attack of -1.5°. This resulted in better agreement on the slat but poorer comparison 
for the main element and flap. Because there was no best USM3D flow solution collection efficiency 
results were generated for both 0° and -1.5° angles-of-attack. 

The collection efficiency comparisons for the three element high-lift model are shown in Figures 22 
to 24. The non-splashing 27 bin 1NS2D and 10 bin USM3D based collection efficiency comparison agree 
well for the angle-of-attack of zero for the slat. The IN S2D non-splashing based calculation generates 
higher collection efficiencies for both the main element and flap than the USM3D based calculation. This 
is probably due to the higher angle-of-attack realized on these elements for the 1NS2D solution than for 
the USM3D solution. It could also be due to the use of the coarser 10 bin distribution used for the 
USM3D based analysis. 

The best agreement for the slat shown in Figure 24 was for the USM3D based flow model at 0° 
angle-of-attack with splashing. The results with splashing for this case resulted in an over-prediction of 
collection efficiency near the lower impingement limit by about 27 percent and an over-prediction of less 
than 5 percent in total collection efficiency. This result is reasonably conservative and much less 
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conservative then the non-splashing results for the 0° angle-of-attack USM3D solution which had an 
over-prediction of total collection efficiency of greater than 1 0 percent. 

Collection efficiency results for the main element are presented in Figure 25. The 1NS2D and 
USM3D based solutions at 0° angle-of-attack over-predict the amount and extent of impingement and the 
number of impingement regions for main element. The experimental data shows two regions of 
impingement on the lower side of the main element. These include a small region near the stagnation 
region and a larger one further back. The analytical solutions reflect impingement in three distinct regions 
one of which is on the upper surface of the main element. It may be that this upper surface impingement 
region was missed during the experiment because strips were not located there. The best agreement 
between experiment and analysis occurred for the -1.5° angle-of-attack based USM3D solution with 
splashing. The extent and magnitude of impingement agree well for this case although the aft 
impingement region occurs further aft for the analytical solution by about 5 cm. 

The multi-element flap collection efficiency results are shown in Figure 26. The analytical results 
under-predicted the amount and extent of impingement and the total collection efficiency for the flap. The 
total collection efficiency was under-predicted by as much as 30 percent for the flap. This could be due to 
an under-prediction in the local angle-of-attack by the flow solver which is evident in the coefficient of 
pressure distribution comparisons for the flap and in the location of the upper impingement limit on the 
flap. The splashing model did not improve the comparisons for flap. A small amount of re-impingement 
on the flap resulted in slightly more mass being added for the USM3D 0° and -1.5° angle-of-attack cases 
but the best agreement was still for the non-splashing INS2D 0° angle-of-attack bases solution. The lack 
of agreement for the flap was probably due to flow-field deficiencies. Past experience has shown that if 
the flow-field matches experiment then the impingement data compares well. Future work needs to be 
done to generate a better flow simulation for this geometry. 

Collection efficiency calculations for the trap wing model are shown in Figures 27 to 31. No 
experimental collection efficiency data was available for this geometry. Calculations for a 10 bin 
distribution were made for a 92 pm MVD drop size indicative of an SLD environment and for a single 
20 pm drop which is more indicative of an Appendix C environment. The 92 pm MVD drop case with 
splashing resulted in about 200 to 500 percent greater total collection efficiency than for the 20 pm drop 
case. The 92 pm MVD drop splashing case yielded about 15 to 30 percent less total collection efficiency 
then the 92 pm MVD drop case without splashing for the trap wing elements shown in Figures 28 to 30. 
The mass loss for the trap wing is similar to that found for the models above. Evident in Figure 3 1 are 
areas near the wing root on the slat and flap which show higher collection efficiencies for the splash 
model due to the re-impingement of the mass splashed from the nose, slat and main elements. 

The DLR-F4 collection efficiency results are presented in Figures 32 to 35. The splashing model 
results in significant mass loss from the nose, windshield and wing. The overall mass loss due to 
splashing for the 92 pm MVD drop case is about 15 to 30 percent for the wing, nose, and windshield. 

A small amount of re-impingement on the base of the windshield from the aircraft nose is evident in 
Figure 33 and 35 (< 5 percent). The 92 pm MVD drop case with splashing yielded about 300 percent 
more mass on the wing and about 800 percent more mass on the nose and windshield than for the 20 pm 
drop case. Areas of higher collection efficiency are evident in Figure 35 for the splashing model near the 
wing root and on the base of the windshield. These are due to re-impingement of mass splashed from the 
nose. 


Conclusion 

The analytical and experimental comparisons for pressure distribution and collection efficiency using 
the NASA Glenn splashing model were in general good. The collection efficiency models in LEWICE 
and LEWICE3D also compared well for the 2D wing models. 

The analytical pressure distributions for all but the multi-element wing compared well with 
experiment. Small differences were noted near the pressure peak and near the trailing edges for most of 


N ASA/TM— 20 11-216945 


9 



the models. These differences were attributable to deficiencies in the inviscid flow-fields and to grid 
resolution. The multi-element wing results proved less then satisfying. The inclusion of viscous effects, 
wind tunnel walls and AOA adjustments did not produce excellent agreement to experimental data. The 
pressure distributions on the slat proved to be the worst while those on the main element and tail were 
reasonable. These differences were enough to cast doubt on the collection efficiency results obtained 
using these flow-fields. It was concluded that more work is require to produce a better multi-element 
flow-field. 

The SLD splashing model generated good agreement with the experimental collection efficiency data 
for all but the multi-element wing. The differences in overall collection efficiency between experiment 
and the analytical splashing model were on the order of about 5 percent for all models except the multi- 
element which was on the order of about 30 percent. The large difference for the multi-element case was 
attributed to flow-field errors. The splashing model fared much better than the standard, non-splashing 
model for which the total collection efficiency differences were as high as 30 percent. 

The collection efficiency results for the trap wing and DLR-F4 cases were similar to those for the 1RT 
test models. The splashing models resulted in a mass loss of about 15 to 30 percent. Some of the splashed 
mass for the cases re-impinged on downstream components. The re-impingement though was in general 
small for these cases and was less than 5 percent of the impinging mass. 
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TABLE 1.— TEN BIN DISTRIBUTION FOR 92 pm 
MVD DROP SIZE GENERATED BY 
WSU SPRAY SYSTEM 


Drop size 
(pm) 

Percent 

LWC 

10.96046 

0.05 

27.59405 

0.1 

51.17964 

0.2 

92.18893 

0.3 

136.2497 

0.2 

179.8884 

0.1 

230.625 

0.03 

83.3301 

0.01 

339.5112 

0.005 

391.7745 

0.005 


TABLE 2.— TWENTY-SEVEN BIN DISTRIBUTION 
FOR 92 pm MVD DROP SIZE GENERATED 
BY WSU SPRAY SYSTEM 


Drop size 
(pm) 

Percent 

LWC 

10.7792 

0.0475 

24.37238 

0.0475 

29.77389 

0.0475 

35.24709 

0.0475 

41.99855 

0.0475 

51.40099 

0.0475 

61.50823 

0.0475 

69.93255 

0.0475 

77.50246 

0.0475 

84.73027 

0.0475 

91.71088 

0.0475 

98.83237 

0.0475 

106.2296 

0.0475 

113.8977 

0.0475 

121.8354 

0.0475 

130.5719 

0.0475 

140.5805 

0.0475 

152.2594 

0.0475 

165.6468 

0.0475 

187.6024 

0.0475 

209.6167 

0.01 

226.1367 

0.01 

247.145 

0.01 

265.6669 

0.005 

290.3102 

0.005 

339.5112 

0.005 

391.7745 

0.005 
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Figure 1. — NASA Glenn Icing Research Tunnel (IRT). Plan view. 



Figure 2. — Installation of spray system used for IRT impingement tests. 
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Figure 3. — Automated reflectometer used to reduce impingement data. 



Figure 4. — Illustration of blotter strip application. 
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(a) IRT installation. 



(b) USM3D surface grid. 
Figure 5. — MS-317 model. 



(a) IRT installation. (b) USM3D surface grid. 

Figure 6.— NLF(1)-0414 model. 
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Figure 7. — Swept 64A008 wing. 



(a) IRT installation. ( c ) USM3D surface 9 rid - 


Figure 8. — Multi-element model. 
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Figure 9. — Surface model for CFD++ trap wing flow model. 



Figure 10. — Surface model for NSU3D DLR-F4 flow model. 
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Figure 1 1 . — Comparison of experimental and analytical coefficient of pressure distribution for MS-317 model. 
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Figure 12. — Comparison of experimental and analytical collection efficiency for MS-317 model. 


N ASA/TM— 20 11-216945 


18 


Coefficienty of Pressure 



-200 -100 0 100 200 300 400 500 600 

Surface Distance From Highlight, mm 


Figure 13. — Comparison of experimental and analytical collection efficiency for MS-317 model. 
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Figure 14. — Comparison of experimental and analytical coefficient of pressure distribution for NLF(1) 0414 model. 
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Figure 15. — Comparison of experimental and analytical collection efficiency for NLF(1) 0414 model. 
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Figure 16. — Comparison of experimental and analytical collection efficiency for NLF(1) 0414 model. 
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Figure 17. — Comparison of experimental and analytical coefficient of pressure for swept NACA 64A008 wing. 



Surface Distance From Highlight, mm 

Figure 18. — Comparison of experimental and analytical collection efficiency for swept NACA 64A008 wing. 
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Figure 19. — Comparison of experimental and analytical collection efficiency for swept NACA 64A008 wing. 



Figure 20. — Multi-element wing cross section. 
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Figure 21. — Comparison of experimental and analytical pressure distribution for multi-element slat. 
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Figure 22. — Comparison of experimental and analytical pressure distribution for multi-element wing main element. 
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Figure 23. — Comparison of experimental and analytical pressure distribution for multi-element model flap. 
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Figure 24. — Comparison of experimental and analytical collection efficiency for multi-element model slat. 
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Figure 25. — Comparison of experimental and analytical collection efficiency for multi-element model main element. 
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Figure 26. — Comparison of experimental and analytical collection efficiency for multi-element wing flap. 
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Figure 27. — Trap wing profiles. 
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Figure 28. — Collection efficiency results for trap wing slat. 
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Figure 29. — Collection efficiency results for trap wing main element. 
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Figure 30. — Collection efficiency results for trap wing flap. 
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Figure 31 . — Collection efficiency results for trap wing. 



Figure 32. — Geometry profiles for DLR-F4 fuselage centerline and wing. 
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Figure 33. — Collection efficiency results for DLR-F4 fuselage centerline. 
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Figure 34. — Collection efficiency results for DLR-F4 wing. 
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Figure 35. — Collection efficiency results for DLR-F4 Airplane. 
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